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MICROWAVE  ENERGY  COUPLING  IN  A  NITROGEN  BREAKDOWN  PLASMA 
1.0  Introduction 

The  pulsed  breakdown  in  air  and  other  gaseous  elements  has  been 
studied  extensively1  with  emphasis  on  the  threshold  power  for  breakdown 
and  its  dependence  on  the  gas  pressure  and  radiation  wavelength.  The 
hydrodynamic  effects  were  first  considered  by  Lin  and  Theofilos.  Their 
calculation  of  the  gas  heating  showed  that  both  high  field  intensity  above 
the  breakdown  intensity  and  increasing  molecular  density  discouraged 
energy  deposition  in  air.  Though  the  preliminary  experimental  results 
confirmed  the  general  features  of  the  simple  theory,  the  observed  pressure 
waves  were  stronger  than  predicted.  Scharfman3  et  al.,  have  observed  that 
a  breakdown  plasma  prevents  the  transmission  of  microwave  enerqy  to  points 
beyond  the  plasma.  Experiments4  recently  performed  at  the  Naval  Research 
Laboratory  (NRL)  showed  that  the  microwaves  are  rapidly  decoupled  from  a 
target  or  plasma  surface.  This  rapid  decoupling  of  the  microwaves  would 
limit  the  amount  of  energy  that  can  be  deposited  into  the  gas.  Computer 
simulations  of  the  NRL  experiments  have  been  performed  at  a  pressure  of  25 
Torr  in  nitrogen.  Throughout  this  work,  we  have  used  the  MINI  code 
developed  at  NRL  to  help  in  understanding  the  basic  interaction  processes 
in  the  experiments. 

2.0  The  Microwave  Nitrogen  Interaction  Code  (MINI) 

The  Mlicrowave  Nitrogen  interaction  code  (MINI)  used  in  our 
studies  is  an  one-dimensional  multi-species,  multi-temperature,  hydro- 
dynamic,  nitrogen  chemistry  and  wave  optics  code.  The  species  followed 
dynamically  were  tfe  ,  N2(X),  N2  (B) ,  N4,  N,  N(2D),  N2(A3Z),  N2(B3it), 

N2(C3it),  and  the  electron  density  (ne).  The  model  calculates  the 
electron  (Te),  vibrational  (Ty)  and  gas  (Tq)  temperatures,  and  is 
space  and  time  dependent  with  microwave  absorption  and  reflection 
considered  in  the  wave  optics  mode.  The  detail  of  the  chemistry  and  the 
wave  optics  aspect  of  the  code  are  described  elsewhere5  and  will  not  be 
repeated  here.  The  electrons  are  described  hydrodynamically  by  the 
equations 
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3.  (n  )  +7* (n  u  )  =  C 
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where  pe  =  ngTe  is  the  electron  pressure,  3g  =  -engug  is  the 
electron  current  density,  Re  =  -meng-)cCu^  -  "u^,)  the 
electron  momentum  transfer  rate,  ue  =  e/mevc  is  the  electron 
mobility,  and  De  =  Te/mgvc  is  the  electron  diffusion  coefficient. 

The  electron  momentum  transfer  frequency  is  vc,  and  me  is  the  mass  of 
the  electron.  The  electron  energy  flux  is  Qe  =  -xeVTg  where  the 
electron  thermal  conductivity  is  <e  =  0ngTe/mevc  and  g  is  a 
constant.  The  terms  on  the  right  hand  side  of  the  energy  equation 
represent  energy  lost  due  to  heat  flow,  work  due  to  momentum  transfer, 
Joule  heating  and  heating  due  to  elastic,  inelastic,  and  chemical 
processes.  The  term  ;e  represents  sources  and  sinks  terms  for  the 
electron  due  to  ioniza'ion,  recombination,  and  attachment.  The  Ee  term 
in  the  electron  energy  equation  includes  energy  losses  due  to  vibrational 
and  electronic  excitation  of  molecular  nitrogen.  The  Ce  and  Eg  terms 
are  discussed  more  fully  in  Ref.  (5).  The  average  velocity  of  molecular 
nitrogen, ”um,  is  related  to  the  diffusion  and  mass-average-velocity 6  of 
the  heavy  particles  by  "1%  =  vm  +  u.  The  diffusion  velocity  is  vm, 
and  u  is  the  mass-average-velocity  of  the  heavy  particles. 
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The  continuity  equation  for  each  of  the  heavy  particle  species 


is 


3t(nl)  +  7*  (n^)  =  C^. 


The  momentum  equation  for  the  heavy  particles  is 

s^ipu)  +  $*{puu  +  n_)  =  p+f  -  i$e 


where  p  =  E^m^n^,  p+  -  Ej^enj,  and  n_  =  tt_  +  p^I  are  the 

mass  density,  positive  charge  density,  and  pressure  tensor.  The  viscous 

portion  of  the  pressure  tensor  is  tt_,  p^  =  TgEjnj  is  the  heavy 

particle  pressure,  and  _I  is  the  unit  tensor.  The  heavy  particle  energy 

equation  is 


3t(eH}  +  ^EHS  +  %  +I*^  =  -5e*f  +  V3e  +  -m^v(Tv)  '  ev(Tq^  +  EH 


£H  =  EiCiniTq  +  “  pu2 


where  i  ‘he  total  energy  density,  c^  is  the 

translation-rotational  specific  heat,  ev(Ty)  =  e{exp(e/Ty)  -  I}"1  is 
the  average  vibrational  energy  per  molecule,7  e  *  0.3  eV  is  the  quantum  of 
vibrational  energy  for  nitrogen,  tv  is  the  characteristic  time  for 
vibrational  relaxation,8  and  E^  is  the  heating  terms  due  to  chemical, 
elastic,  and  inelastic  processes.  Again  a  more  detailed  discussion  of 
can  be  found  in  Ref.  (5).  The  heavy  particle  energy  flux  is  QH  = 
anc*  *h  thermal  conductivity  of  the  gas,  including 

translational  and  rotational  contribution.8  The  average  and  diffusion 
velocity  of  species  i  are  related  byH^  =  vj  +  IT.  The  diffusion 
velocity  of  the  ith  ion  species  is  "v\  and  TT  is  the  mass-average-velocity 
of  the  heavy  particles 


u  = 


imiVi  51  V 
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The  vibrational  energy  equation  may  be  written  approximately  as 


3  (n  e  (T  )}  +  $*{n  e  ( T  ) u  +  $  }  =  - —  {e  (T  )  -  e  (T  )}  +  E 

tlmvv)  lmvvm  v1  t  lvv  vgJ  v 


where  Qy  =  -tcvVTv  is  the  vibrational  energy  flux  and  <v  is  the 
vibrational  conductivity.  The  vibrational  temperature  can  be  calculated 
assuming  that  the  nitrogen  molecule  Is  a  harmonic  oscillator  and  that  the 
vibrational  levels  have  a  Boltzman  distribution.  Hence,  the  vibrational 
energy  source  term  Ev  is 


E 

v 


en  n  f  1 
m  el 


exp(-e/Tv)}y®=1vXv{l 


exp[ev(Tv  -  T^/TjJ} 


where  Xy  is  the  excitation  rate  coefficient10  for  the  vth  vibrational 
level  obtained  from  the  experimentally  measured  cross-sections.11  The 
vibrational  energy  spacing  is  e  =  0.3  eV  for  nitrogen. 


The  equations  solved  in  MINI  are: 

3.  fnl+3  (n  u  )  i  C 
t  ^  eJ  x  ^  e  e'  e 
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(2) 
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(6) 


(7) 


3 .  { n  e  ( T  1}  +  3  {nc  f  T  )  u }  =  -  —  { e  ( T  )  -e  ( T  )}  + 
m  v1-  v'1  xl  m  v^  v;  1  t  1  v'-  v'  v^  q'J 


where  the  Kinematic  and  bulk  viscosity  coefficients  are  y  and  y^.  The 
transport  coefficients  for  the  viscosity  and  thermal  conductivity  are 
taken  from  Bird.12 

The  dc  portion  of  the  electric  field  and  the  diffusion  veloci¬ 
ties  are  determined  by  demanding  the  total  current 

\  =  £  *  y.3.  =  0  (8) 

T  e  -  j.  i 

where  0^  =  Z^en^u^  is  the  current  density  of  the  ith  ion  species. 

In  steady  state,  each  of  the  ion  species  satisfies  an  equation 

u.  =  u  +  y,?  -  D.V(lnn.)  -  D  V(lnT  )  (9) 

I  m  i  i  l  i  g 

where  the  ion  diffusion  coefficient  and  mobility  are  related  to  the  colli¬ 
sion  frequency  vj  by  Dj  =  T  /m^v^  and  y^  =  Z^e/m^v^. 

The  mass  of  the  ith  ion  species  is  mj.  Summing  Eq.  (9)  for  all  ion 
species  and  using  Eq.  (9)  gives  an  equation  for  the  ambipolar  electric 
field.  Once  the  ambipolar  field  is  found,  Eqs.  (9)  and  (2)  are  used  to 
solve  for  the  flow  velocities  of  the  ions  and  the  electron. 

The  Eulerian  difference  equations  of  Eqs.  (1)  -  (7)  are  solved 

i  3 

using  the  algorithm  of  Rubin  and  Burstein.  The  scheme  is  simple,  second 
order  in  accuracy,  and  is  relatively  stable  while  centering  the  dis¬ 
sipation  terms.  The  hydrodynamic  and  chemistry  portions  are  calculated 
separately  in  a  cycle,  assuming  the  flow  velocities  of  the  fluid  are 
weak.  The  hydrodynamics  is  determined  from  the  updated  values  of  the 
chemistry.  Hence,  the  hydrodynamics  is  considered  as  a  small  correction 
to  the  chemistry. 
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3.0  Cascade  Ionization  of  Nitrogen  by  a  Microwave  Pulse 

Microwave  gas  breakdown  begins  when  a  small  number  of  priming 
electrons  acguire  sufficient  energy  from  the  electric  field  to  ionize  the 
gas.  These  priming  electrons  gain  energy  when  the  oscillatory  energy  of 
the  electrons  in  the  electric  field  is  randomized  by  collisions  into 
thermal  energy.  Once  an  electron  has  an  energy  in  excess  of  the  ioni¬ 
zation  energy  of  the  gas,  the  electron  can  easily  ionize  the  gas,  result¬ 
ing  in  the  generation  of  two  low  energy  electrons.  The  electron  avalanche 
can  be  described  by 


ng(t)  =  ne(0)exp(t/xB) 


(10) 


where  ne(t)  Is  the  electron  density  at  time  t,  and  Tg  is  the 
breakdown  time.  Physically,  tg  is  the  time  required  for  an  electron 
to  gain  the  energy  needed  to  ionize  the  gas.14  Fig.  1  shows  the  breakdown 
time,  Tg(sec),  and  the  electron  temperature,  Te(eV),  for  conditions 
typical  of  microwave  breakdown  of  N2 .  A  high  radiation  field  would  pro¬ 
mote  the  electron  avalanche  by  increasing  the  enerqy  imparted  to  the 
electron  per  collision  with  the  molecule.  A  high  molecular  density,  how¬ 
ever,  requires  higher  intensities  for  the  avalance  breakdown.  \t  moderate 
pressures  above  25  Torr,  the  dominant  energy  loss  mechanism  for  the 
electrons  is  the  excitation  of  the  electronic  and  vibrational  modes  of  the 
gas.  The  root-mean-square  electric  field  to  density  ratio,  E/nm 
(V-cm  ),  does  not  uniquely  define  the  absorbed  energy  for  an  ac  electric 
field.  The  striped  region  in  Fig.  1  is  bounded  by  the  constant  intensity 
line  (I)  and  the  constant  pressure  line  (P).  The  rapid  variation  in 
Tg  with  increasing  Intensity  allows  a  breakdown  threshold,  Ig,  to 
be  assigned  at  the  onset  of  breakdown.  The  ability  to  couple  microwave 
energy  to  a  breakdown  gas  depends  sensitively  on  the  breakdown  time, 
tr,  and  the  electron  density. 
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(sec) 


The  efficiency  of  converting  microwave  energy  to  the  trans¬ 
lational  energy  of  the  gas  can  be  enhanced  by  coupling  the  microwaves  to  a 
preformed  plasma.  Plasma  maintenance  of  a  preformed  plasma  by  microwaves 
ailows  the  microwave  deposition  to  be  localized  while  maintaining  a  high 
electron  density.  In  short  pulse  breakdown  experiments  with  pulse 
lengths,  tp  <  1  ysec,  a  high  electron  density  is  a  necessary  condition 
for  rapid  heating  of  the  gas.  Plasma  maintenance  requires  the  breakdown 
time  T0  >  when  n„(0)  >  1.  If  the  local  electric  field  is  not 
carefully  matched  to  the  instantaneous  plasma  conditions,  the  gas  heating 
will  be  self-iimited  by  the  formation  of  a  non-hydrodynamic  ionization 
front.  Gas  breakdown  off  a  reflecting  surface  without  a  preformed  plasma 
will  not  result  in  a  local  deposition  of  the  microwave  energy.  Since  the 
microwave  deposition  is  determined  self-consistently  with  the  electron 
density  profile,  the  resulting  plasma  density  cannot  rapidly  heat  the  gas 
in  the  short  pulse  experiments. 

3.1  Gas  Breakdown  Near  a  Reflecting  Surface 

A  plane  wave  incident  normal  to  a  reflecting  surface  would  re¬ 
quire  less  power  for  gas  breakdown.  The  constructive  interference  of  the 
incident  and  reflected  microwave  results  in  the  characteristic  standing 
wave  pattern  for  the  electric  field  and  a  factor  of  four  increase  In  the 
field  intensity.  If  the  intensity  is  above  the  breakdown  threshold, 
breakdown  of  the  gas  occurs  instantaneously  at  the  peak  nodes  of  the  elec¬ 
tric  field.  As  the  electron  density  rises,  these  "pancakes"  of  electrons 
begin  to  attenuate  the  radiation.  Pancakes  closest  to  the  reflecting  sur¬ 
face  rise  more  rapidly  early  in  the  pulse  (Fig.  2).  As  the  electron 
density  continues  to  rise,  attenuation  of  the  microwave  causes  the  pan¬ 
cakes  closest  to  the  surface  to  fall  in  density  below  that  of  pancakes  far 
from  the  surface.  Shorter  systems  with  fewer  pancakes  decouple  more 
rapidly  than  longer  systems.  However,  the  total  absorbed  energy  remains 
essentially  constant,  independent  of  the  system  size  for  a  fixed  pressure 
and  pulse  length.15 
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Fig.  2  —  Gas  breakdown  near  a  reflecting  surface  showing 
the  electron  density  pancakes  for  (a)  one-and-a  half,  s..i 
(b)  three  wavelengths  planar  system.  The  critical  density  is 
nr  =  1.5  X  1013  cm*3  and  k  =  7.33  cm'1. 


The  finite  length  simulations  represent  the  planar  region  near 
the  surface  of  a  focused  microwave  system  where  the  distance  from  the  sur¬ 
face  is  less  than  the  radius  of  the  spotsize.  Fig.  2  shows  two  similar 
simulations  with  I  =  6.25  kW/cm2  for  an  one-and-a-half  and  three  wave¬ 
length  system  without  hydrodynamic  effects.  The  plasma  absorbs 
approximately  20%  of  the  total  microwave  energy  with  an  asymptotic 
absorption  efficiency  of  80%.  The  absorption  is  regulated  by  the  electron 
density.  Shorter  systems  have  a  much  higher  electron  density  than  larger 
systems  and  can  more  readily  heat  the  gas  in  the  short  pulse  experiments. 
The  gas  is  not  expected  to  be  heated  more  than  the  observed  4%  above  the 
ambient  temperature  in  a  weakly  focused  system.  The  absorbed  energy  is 
distributed  over  a  larger  number  of  pancakes  with  a  very  low  electron  den¬ 
sity.  As  will  be  discussed  later,  the  low  electron  density  of  Fig.  2 
results  in  a  long  heating  time  for  the  gas  of  =  6.5  ysec  where 
h  Tg( dTg/dt) -1 .  The  time  required  to  heat  the  gas  to  an 
appreciable  temperature  in  a  weakly  focused  system  would  require  pulse 
lengths  of  Tp  >  6.5  ysec.  However,  even  for  long  pulses,  a  significant 
heating  of  the  gas  is  never  realized.  Later  in  the  pulse,  the  pancakes 
furthest  away  from  the  reflecting  surface  continue  to  rise  in  density. 
Ultimately,  a  single  pancake  will  have  sufficient  density  to  reflect  the 
radiation.  This  pankcake  acts  as  a  secondary  reflecting  surface  and 
decouples  the  radiation  from  the  primary  metallic  surface.  Gas  breakdown 
experiments  using  a  reflecting  surface  show  this  behavior.  Similarly, 
free  nitrogen  breakdown  simulations  without  a  surface  using  a  focused 
system  also  show  that  the  radiation  is  decoupled  from  the  initial 
breakdown  site. 


3.2  Maintenance  of  a  Preformed  Plasma  by  Microwaves 

Microwave  coupling  to  a  pre- ionized  plasma  near  a  metallic  sur¬ 
face  show  some  general  character istics.  Though  the  plasma  absorbs  the 
microwave  energy,  the  absorption  cannot  be  localized.  A  traveling  ioni¬ 
zation  front  is  observed  to  proceed  toward  the  radiation  source.  Fig.  3 
shows  the  time  history  of  the  ionization  front  in  the  strong 
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ELECTRON  DENSITY  (ne/nc) 


POSITION  (ky) 

Fig.  3  —  Time  history  of  the  ionization  front  in  a  planar 
system  showing  the  wave  motion  toward  the  radiation 
source.  The  microwaves  are  incident  from  the  left  onto 
a  pre-ionized  plasma  near  a  reflecting  surface  on  the  right 
boundary.  The  normalization  parameters  are  the  same  as 
in  Fig.  2. 
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mismatch  case  I  >  Ig  and  Xp  >  x^.  The  critical  time  x^  is  the 
time  required  for  the  electron  avalanche  to  reach  the  critical  density  nc 
=  (o2me/4ire2  =  1013/X2(cm)  where  A  is  the  wavelength  of  the  incident 
radiation.  In  this  simulation  I  =  12.5  kW/cm2  >  Ig  =  3  kW/cm2,  Xp  =  0.5 
usee,  x^,  =  0.16  usee,  and  all  hydrodynamic  effects  are  suppresed.  The 
initial  electron  density  (pre-ionized  plasma)  was  initialized  with  a 
symmetric  profile  with  a  peak  density  of  0.1nc  and  FWKW  of  0. lx.  In  Fig. 

3,  the  peak  density  is  determined  by  Eq.  (10).  However,  the  breakdown  time, 
x D ,  is  a  complicated  function  of  the  incident  and  reflected  microwaves. 

It  is  clear  that  the  small  reflectivity  (R  <  10%)  is  important  in 

determining  the  propagation  speed  since  as  Fig.  3  shows,  the  separation 

distance  between  peaks  is  A/4.  The  ionization  front  proceeds  toward  the 
source  at  a  characteristic  speed  D(t)  =  A/4Xp(t).  The  maximum  density 
is  decreasing  with  time  while  the  total  absorption  continues  to  increase 
with  time.  The  asymptotic  absorption  efficiency  is  approximately  95%.  The 
movement  of  the  ionization  front  limits  the  time  the  gas-plasma  systems 
samples  the  radiation  field.  The  sampling  time,  x^,  Is  determined  by 

a(cm)  =  0(t)dt  *  0{0)x<.  where  D(0)  =  A/4Xp(0)  a  1.6  x  106  cm/sec  is 

the  early  breakdown  speed.  Estimating  the  absorption  length  as  a(cm)  =  A/2tt 
gives  a  sampling  time  of  x  g  *  x^.  In  the  pressure  regime  P  >  25  Torr, 
the  gas  is  primarily  heated  through  the  quenching  of  the  electronic  states  by 
molecular  and  atomic  nitrogen.  The  fractional  change  in  the  gas  temperature 
is  £g  =  (ATg/Tg)Max  =  xsEaBQaBnans£a/cvnmTg  where 

(^0  is  the  quenching  rate  of  the  N2  triplet  state  na  by  the  quenchant 
ng .  The  specific  heat  of  nitrogen  is  cy  and  Ea  is  the  energy  level  of 
the  state.  As  an  example,  the  N2(C  it)  state  is  quenched  by  N2  with  a  rate 
coefficient16  of  1.2  x  10"11  cm3 /sec.  The  fractional  change  in  the  gas 
temperature,  Eg»  is  approximately  21%  from  the  ambient  temperature. 

The  region  of  the  maximum  gas  heating  is  at  the  Initial  breakdown 
site  (  X/4  from  the  surface).  The  electric  field  is  less  than 


1%  of  the  incident  time  average  electric  field  at  the  maximum  gas  tempera¬ 
ture  (Fig.  4a).  The  plasma  absorbs  approximately  70%  of  the  Incident 
microwave  energy  and  has  an  electron  temperature  of  less  than  2.5  eV 
(Fig.  4b).  Maximum  gas  heating  requires  the  sampling  time  and  the  elec¬ 
tron  temperature  to  be  as  large  as  possible  without  decoupling  the 
system.  A  high  electron  temperature  populates  the  electronic  states  which 
are  rapidly  quenched  to  heat  the  gas.  Approximately  4%  and  40%  of  the  in¬ 
cident  energy  goes  into  heating  the  gas  and  exciting  the  vibrational 
states.  A  similar  simulation  with  the  peak  power  reduced  to  I  =  6.25  kW/ 
cm  shows  no  movement  of  the  initial  ionization  region  while  increasing 
the  fractional  change  in  the  gas  temperature  to  27%.  Further  irradiance 
would  increase  the  gas  temperature.  However,  these  results  are  sensitive 
to  the  hydrodynamic  response  of  the  plasma.  Again,  a  similar  simulation 
with  hydrodynamic  effects  showed  a  reduction  in  the  gas  and  vibrational 
temperature.  The  microwave  plasma  system  again  decoupled;  however,  the 
hydrodynamic  motion  reduced  the  energy  input  into  the  gas  and  vibrational 
energy  to  0.1%  and  26%  respectively.  The  electron  density  scale  lengths 
are  sharp,  L/X  >  10"2  where  L  =  ne(dne/dx)_1 .  The  hydrodynamic  motion 
of  the  electrons  tends  to  reduce  the  electron  temperature  and  density  in 
regions  of  high  density.  The  slow  gas  flow  (<  300  cm/sec)  further  reduces 
the  electron  temperature  by  convecting  the  vibrational  energy  away  from 
regions  of  high  electron  density.  The  small  temperature  and  density 
changes  translate  to  greater  than  10%  reduction  in  the  population  of  the 
electronic  states.  Since  quenching  of  the  electronic  state  is  the  domi¬ 
nant  mechanism  for  rapidly  heating  the  gas,  the  gas  heating  is  negligibly 
small.  It  Is  clear  that  the  motion  of  the  ionization  front  will  certainly 
limit  the  gas  heating. 

3.3  Gas  Breakdown  in  a  Focused  Microwave  System 

Simulations  of  a  focused  microwave  system  were  also  performed. 

If  the  field  intensity  in  a  focused  system  greatly  exceeds  the  breakdown 
threshold  at  the  focal  spot,  a  breakdown  wave  propagates  toward  the 
source.17  The  motion  of  the  microwave  plasma  absorbing  layer  is  not  a 
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E  SQUARE  (ARB) 


(a)  (b) 

Fig.  4  —  (a)  The  time  averaged  electric  field  square  and 
(b)  the  electron  temperature  vs  position.  The  normalization 
parameters  are:  Te0  =  0.25  eV,  k  =  7.33  cm'1,  and  the 
arbitrary  units  of  <E2>  =  1  corresponds  to  an  average 
intensity  of  9.4  kW/cm2. 


simple  breakdown  wave,  since  the  reflected  microwave  is  Important  in 
determining  the  propagation  speed.  Fig.  5  shows  the  time  history  of  the 
electron  density  profile  during  the  course  of  a  tp  =  0.5  psec  pulse  in  a 
focused  system.  The  spherical  wave  equation  is  solved  for  the  electric 
field  to  model  the  convergence  of  the  field.  The  peak  incident  field  is  I 
=  10  kW/cm2  with  a  system  aspect  ratio  of  2.5.  The  average  power  at  focus 
is  31.25  kW/cm  ,  and  all  hydrodynamic  effects  are  suppressed.  The  very 
rapid  early  time  speed  of  the  absorbing  wave,  D  *  5  x  106  cm/sec,  gives  a 
sampling  time  of  Xg  =  1.7  x  10"  sec.  However,  the  gas  continues  to 
be  heated  after  the  passage  of  the  field  by  the  de-excitation  of  the  elec¬ 
tronic  states.  The  maximum  gas  and  vibrational  temperatures  are  eg  * 

10%  and  ey  =  260%  respectively.  The  quarter-wavelength  fine  structure 
in  the  density  profile  shows  the  effect  of  the  reflected  light.  The  time 
average  reflection  and  absorption  coefficient  are  19%  and  68%  respect¬ 
ively.  A  snapshot  of  the  electric  field  and  electron  temperature  (Fig.  6) 
shows  the  lack  of  field  penetration  (a  -  0.1X)  and  the  cold  plasma 
temperature,  Te  <  3.3  eV.  The  gas  and  vibrational  temperatures  again 
are  reduced  when  hydrodynamic  effects  are  included  in  the  simulations. 

The  gas  is  expected  to  be  heated,  at  most,  to  ?0%  of  the  initial  value 
during  the  course  of  a  0.5  psec  pulse. 


The  propagation  velocity  of  the  breakdown  wave  in  a  focused 
microwave  system  can  be  calculated.  Simulations  show  the  absorption  is 
localized  at  the  head  of  the  breakdown  wave.  The  reflection  coefficient 
is  relatively  constant  after  the  Initial  breakdown  (t  >  x^).  The 
breakdown  time,  x^,  is  inversely  proportional  to  the  light  flux, 

1/tg  <*  (1  +  2Re^  +  £2)I»  where  the  reflection  coefficient  is  c} .  The 
intensity  in  a  conically  focused  system  can  be  written  as 


I(X,t)  « 


(1  +  2Rec  +  e2) 


x  $  (t ) 


where  $(t)  is  a  form  function  characterizing  the  pulse.  The  radius  of  the 
radiation  channel  at  a  position  X  on  the  channel  axis  is  r  =  r0  +  Xtane. 
The  minumum  radius  of  the  waist  at  the  focal  point  is  r0  and  9  is  the 
half-angle  of  the  radiation  column,  following  Raizer,17  the  propagation 
velocity  of  the  breakdown  wave  is  approximately  D  =  r0/x^,tan9.  Under 
the  assumptions  made,  the  velocity  D  would  be  an  upper  bound  on  the 
propagation  velocity  of  the  breakdown  wave. 
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Fig.  6  —  (a)  The  time  averaged  electric  field  square  and 
(b)  the  electron  temperature  vs  position  for  the  simulation  of 
Fig.  5.  The  normalization  parameters  are  the  same  as  in  Fig.  4. 
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4.0  Simulation  of  the  NRL  Experiments 

We  compare  the  simulation  results  with  the  NRL  experiment.  The 
NRL  35  GHz  112  kW  focused  gyrotron  system  gives  an  average  power  at  the 
focal  spot  of  33  kW/cm  .  The  polarization  of  the  microwave  is  such  that 
the  electric  field  is  perpendicular  to  the  plasma  gradient  (S  polari¬ 
zation).  The  average  aspect  ratio  of  the  experiment  qives  an  equivalent 
simulation  aspect  ratio  of  2.5.  The  free  nitrogen  breakdown  experiments 
were  performed  at  25  Torr  with  a  100  pps  repetition  rate  at  tp  =  1.5 
nsec.  The  minimum  radius  of  the  waist  at  the  focal  spot  is  r0  =  0.75  cm 
and  the  half-angle  of  the  radiation  column  satisfies  t  n9  =  0.35.  Experi¬ 
mentally,  framing  camera  data  (Fig.  7)  gives  a  propagation  velocity  of  D  » 
4.6  x  106  cm/sec.  The  width  of  the  absorbing  wave  is  approximately  1.2 
cm.  The  measured  absolute  intensity  of  the  second  positive  band,  N2  2P 

(0-0)  at  3371  A  is  1.4  x  1021  eV/cm2-sec.  The  calculated  spatial  width  of 

2  1 

the  3371  A  light  emission  is  0.8  cm  with  a  peak  intensity  of  1.8  x  10 
eV/cm  -sec.  As  mentioned  previously,  the  calculated  speed  of  the  break¬ 
down  wave  is  5  x  106  cm/sec.  Experimentally ,  the  inferred  breakdown  time 
is  tb  =  2.5  x  TO'8  sec  or  =  7.4  x  10'7  sec.  This  is  to  be 
compared  to  the  calculated  breakdown  time  early  in  the  pulse,  tg  »  4.5 
x  10"9  sec,  and  late  in  the  pulse,  ig  a  6  x  10'8  sec.  The  simulation 
results  are  felt  to  be  in  reasonable  agreement  with  the  experiment. 

5.0  Conclusion 

In  conclusion,  computer  simulation  of  microwave  coupling  to  a 
weakly  ionized  nitrogen  qas  plasma  shows  excellent  absorption  of  the 
microwave  energy.  The  electronic  states  of  nitrogen  excited  by  electron 
impact  are  rapidly  quenched  by  molecular  and  atomic  nitrogen  to  heat  the 
gas  at  pressures  above  25  Torr.  However,  high  power  irradiance  can  easily 
produce  a  non-hydrodynamic  ionization  front  which  can  severely  limit  the 
heating  of  the  gas.  Microwave  qas  breakdown  experiments  with  pulse 
lengths  less  than  one  microsecond  are  expected  to  show  only  a  slight  heat¬ 
ing  of  the  gas. 
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Fig.  7  —  Framing  namera  data  showing  the  motion  of  the 
ionization  front  in  a  focused  free  nitrogen  breakdown  experi¬ 
ment.  The  top  and  bottom  photographs  were  respectively 
taken  0.2-  and  0.1-  //sec  after  the  first  observable  break¬ 
down. 
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